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Abstract 

We investigate the problem of computing tensor product multiplicities for 
complex semisimple Lie algebras. Even though computing these num- 
bers is #P-hard in general, we show that if the rank of the Lie algebra is 
assumed fixed, then there is a polynomial time algorithm, based on count- 
ing the lattice points in polytopes. In fact, for Lie algebras of type Ar, 
there is an algorithm, based on the ellipsoid algorithm, to decide when 
the coefficients are nonzero in polynomial time for arbitrary rank. Our 
experiments show that the lattice point algorithm is superior in practice 
to the standard techniques for computing multiplicities when the weights 
have large entries but small rank. Using an implementation of this algo- 
rithm, we provide experimental evidence for conjectured generalizations 
of the saturation property of Littlewood-Richardson coefficients. One of 
these conjectures seems to be valid for types _B„, C„, and _D„. 

1 Introduction 



Given highest weights A, /x, and v for a finite dimensional complex semisimple 
Lie algebra, we denote by C^^ the multiplicity of the irreducible representation 
Vi, in the tensor product of V\ and V^; that is, we write 

^a®^m = 0Ca^K. (1) 

V 

In general, the numbers are called Clebsch-Gordan coefficients. In the 
specific case of type A,. Lie algebras, the values C^^ defined in equation ^ are 
called Littlewood-Richardson coefficients. When we are specifically discussing 
the type A^ case, we will adhere to convention and write c^^ for C^^. 
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first author is also grateful for support received from the Alexander von Humboldt foundation 
and a UC Davis Chancellor's fellow award. 
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The concrete computation of Clebsch-Gordan coefficients (sometimes known 
as the Clebsch-Gordan problem \21\ ) has attracted a lot of attention from not 
only representation theorists, but also from physicists, who employ them in the 
study of quantum mechanics (e.g. |lU[l37j '). The importance of these coefficients 
is also evidenced by their widespread appearance in other fields of mathematics 
besides representation theory. For example, the Littlewood-Richardson coef- 
ficients appear in combinatorics via symmetric functions and in enumerative 
algebraic geometry via Schubert varieties and Grassmannians (see for instance 
|3H1E21). More recently, Clebsch-Gordan coefficients are playing an important 
role on the study of P vs. NP (see UHl)- Very recently, Narayanan has proved 
that the computations of Clebsch-Gordan coefficients is in general #P-complete 
123. Here are our contributions: 

(1) We combine the lattice point enumeration algorithm of Barvinok j3] with 
polyhedral results due to Knutson and Tao and Berenstein and Zelevinsky 
|S] in the polyhedral realization of Clebsch-Gordan coefficients to produce a 
new algorithm for computing these coefficients. Our main theoretical result is: 

Theorem 1.1. For fixed rankr, if q is a complex semisimple Lie algebra of rank 
r, then one can compute Clebsch-Gordan coefficients for g in time polynomial 
in the input size of the defining weights. 

Moreover, as a consequence of the polynomiality of linear programming and 
the saturation property of Lie Algebras of type A,., deciding whether c"^^ = can 
be done in polynomial time, even when the rank is not fixed. 

(2) We implemented the algorithm for types Ar, Br, Cr, and Dr (the so-called 
"classical" Lie Algebras) using the software packages LattE and Maple. In 
many instances, our implementation performs faster than standard methods, 
such as those implemented in the software LiE. Our software is freely available 
at |http : //math . ucdavis . edu/~tmcal| 

(3) Via computer experiments, we explored general properties satisfied by the 
Clebsch-Gordan coefficients for the classical Lie algebras under the operation of 
stretching of multiplicities in the sense of |18| . On the basis of abundant exper- 
imental evidence, we propose two conjectured generalizations of the Saturation 
Theorem of Knutson and Tao One of them, which applies to all of the 
classical root systems, is an extension of earlier work by King et al. ([ISp. 

Organization of the paper: In Section|21 after a review of the background mate- 
rial, we prove Theorem ll.il Section |31 explains our experiments comparing our 
software, a mixture of Maple and LattE |13| . with LiE. In Sectional we present 
the two conjectures, both of which, if true, would generalize the Saturation 
Theorem of Knutson and Tao. 
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2 Clebsch— Gordan coefficients: Polyhedral Al- 
gorithms 



As stated in the introduction, we are interested in the problem of efficiently 
computing C^^ in the tensor product expansion V\ = ©i^^v^'^' 
appears that the most common method used to compute the Clebsch-Gordan 
coefficients is based on Klimyk's formula (see lemma below). For example, it is 
used in LiE [3^1 and the Maple [2^] packages Coxeter and Weyl |34| . 

Lemma 2.1. fllGf . Exercise. 24-9) Fix a complex semisimple Lie algebra g, 
and let W be the associated Weyl group. For each weight v of Q, let sgn^v) 
denote the parity of the minimum length of an element cr S 22J such that (7(1^) is 
a highest weight, and let {v} denote that highest weight. Let 6 be one-half the 
sum of the positive simple roots of q. Finally, for each highest weight A of g, let 
K\i, be the multiplicity of v inV\. 

Then, given highest weights A and fi of g, we have that 

14 = Kxe sgn(£ + 1^1 + 5)V[,+^+s^,_s, 

s 

where the sum is over weights s of g with trivial stabilizer subgroup in 20. 

Implementations of Klimyk's algorithm begin by computing the weight spaces 
appearing with nonzero multiplicity in the representation Vx. Then, for each 
such weight e with trivial stabilizer, one computes the Weyl group orbit of 
e + n + S. One then finds the dominant member of the orbit and notes the 
number I of reflections needed to reach it. Finally, one adds (— l)'i^A6 to the 
multiplicity of V{^+i_,+s}-5- 

From the point of view of computational complexity, Klimyk's algorithm has 
two main disadvantages. First, it requires the computation of weight space mul- 
tiplicities, which is in general a #P-hard problem [27]. The second disadvantage 
is that the algorithm requires visiting all of the orbit members, which can be an 
exponentially large set. Indeed, Klimyk's formula above is exponential in the 
size of the input weights. Thus, in practice, the sizes of A, ^, and 1/ usually need 
to be small. One can then ask for an algorithm that behaves well as the sizes of 
the input weights increase, at least if some other parameter is fixed. Stembridge 
also raised the challenge of crafting algorithms based on geometric ideas such as 
Littelmann's paths or Kashiwara's crystal bases |23l 1^ (see comment on page 
29, section 7, of As we see below, there is such an algorithm, based on 

the polyhedral geometry of the Clebsch-Gordan coefficients. 

In 1992, Berenstein and Zelevinsky presented a combinatorial interpreta- 
tion of the Littlewood-Riehardson coefficients as the number of lattice points 
in members of a certain family of polytopes |5]. In 1998, Knutson and Tao 
introduced another family, the hive polytopes, which they used to prove the 
Saturation Theorem (see jj^l and Section^. Each of the polytopes presented 
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by Berenstein and Zelevinsky in 1992 is the image under an injective lattice- 
preserving linear map of a hive polytope |28| . Therefore d^^ equals the number 
of integer lattice points in a corresponding hive polytope H^^. Finally, in 2001, 
Berenstein and Zelevinsky [H] introduced polytopes which enumerate Clebsch- 
Gordan coefficients for any finite dimensional complex semisimple Lie algebra. 
We refer to this last family of polytopes as the BZ~polytopes. We now give the 
exact definition of the hive polytopes. These polytopes exist in the polyhedral 
cone of hive patterns, which we now define. 

Definition 2.2. Fix r e Z>o and let H = k) G Z^.^ : i + j + fc = r}. A 

hive pattern is a map 

/i:7Y^M>o, {i,j,k)^hijk, 
satisfying the rhombus inequalities: 

• /ii,i-l,fc+l + ^i-lj" + l,fc < hiji; + /li-lj-,fc+l, 

• /ii+lj'-l,fc + < hijk + 

for {i,j, k) e Ti, i,j > 1. 

Equivalently, a hive pattern of rank r is a triangular array of real numbers 

^1,0, r-l ^0,1, r-l 

^2,0,r-2 /ll,l,r-2 /l0,2,r-2 



such that, in every "little rhombus" of entries 



c b a c 

a b 

ad d b 

d 



we have a + b > c + d. 
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Here is example of a hive pattern with r = 5: 



5 8 
8 12 13 

11 15 17 18 

12 16 18 20 20 

Recah that when g is of type Ar, so that g ^ 51t.+i(C) for some r > 2, the 
highest weights are (with respect to the canonical basis) partitions of length r, 
i.e., sequences A of integers Ai > • • ■ > A„ > 0. Wc write |A| for J2i the size 
of the partition A. 

Definition 2.3. Given partitions X,ii,i^ G Z>oi the hive polytope H'(^ is the 
set of hive patterns with boundary entries fixed as below. 



i^i = • • = Ai 

Vi + 1/2 = * • • = Ai + A2 

+ 1/2 + 1^3 = • • • • = Ai + A2 + A3 



• • • = |A| 

^ ^ ^ 

X X X 

V y V 

X X 



Note that, for fixed r, the input size of a hive polytope grows linearly 
with the input sizes of the weights A, /i, and v. We will need the following result: 

Lemma 2.4. fl2(Jf ) The Littlewood-Richardson coefficient c\^ equals the num- 
ber of integer lattice points in the hive polytope H^^ . 

Unfortunately, the description of the BZ-polytopes is more involved than 
that of the hive polytopes above. Therefore, we refer the reader to Theorems 2.5 
and 2.6 of [S], which give their description as systems of linear equalities and in- 
equalities in terms of the root systems Br, Cr, and Dr- The reader may also view 
the contents of our maple notebook, available at'http: //math.ucdavis . edu/~tmcal^ 
for completely explicit descriptions of the necessary inequalities. The specific 
properties of the BZ-polytopes that we need to prove our theorem are (1) for 
fixed rank r, the dimensions of the BZ-polytopes are bounded above by a con- 
stant, (2) the input size of a BZ-polytope grows linearly with the input sizes of 
the weights A, /i, and i/, and (3) the following result describing the relationship 
between BZ-polytopes and Clebsch-Gordan coefficients: 
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Lemma 2.5. (Theorem 2.4 of J^) Fix a finite dimensional complex semi-simple 
Lie algebra g and a triple of highest weights {X,^,v) for q. Then the Clebsch- 
Gordan coefficient C^^ equals the number of integer lattice points in the corre- 
sponding BZ-polytope. 

Finally, the last ingredient necessary is A. Barvinok's algorithm for counting 
lattice points in polytopes in polynomial time for fixed dimension. Several de- 
tailed descriptions of the algorithm in Lemma l2. Bl are now available the literature 
(see a-nd all the references therein). 

Lemma 2.6. fWlj) Fix d e Z>o. Then, given a system of equalities and inequal- 
ities defining a rational convex polytope P C K*^, we can compute ^{P f) Z'^) in 
time polynomial in the input size of the polytope. 

Proof of Theorem \l.l\ First, if we fix the rank of the Lie algebra, then we fix an 
upper bound on the dimension of the hive or BZ polytope. Moreover, the input 
sizes of these polytopes grow linearly with the input sizes of the weights. Thus, 
by Barvinok's theorem (Lemma 12.61 stated above) , their lattice points can be 
computed in time polynomial in the input sizes of the weights. Therefore, the 
theorem follows by Lemmas I2.4l and l2.5l 

For the second part of the theorem regarding type Ar, the hive polytopes 
provide a very fast method for determining whether c\^ — 0. According to the 
Saturation Theorem (see Section^)), c\^ = if and only if the corresponding hive 
polytope is empty. Hence, it suffices to check whether the system of inequalities 
defining the hive polytope is feasible, which can be done in polynomial time for 
arbitrary dimension as a corollary of the polynomiality of linear programming 
via Khachian's ellipsoid algorithm (see |3U| ) . 

It is useful to notice that every Kostka number if is a Littlewood-Richardson 
coefficient for some choice of highest weights. One short bijection is given by 



For those familiar with the enumeration of semi-standard Young tableaux 
and Littlcwood-Richardson tableaux by Kostka numbers and Littlcwood-Rich- 
ardson coefficients respectively (see, e.g., (32|), the bijection is straightforward: 
Given a semi-standard Young tableau Y with shape A and content /i, construct 
a Littlewood-Richardson tableau L with shape r/tr and content A by filling the 
boxes as follows. Start with a skew Young diagram D with shape t/ct. For 
J = 1, . . . , r, and for each entry i in the jth row of Y , place a j in the ith row 
of D. Let L be the tableau produced by filling the boxes of D in this fashion. 
(See Figure^) It is not hard to see that, under this map, the column-strictness 
condition on Y is equivalent to the lattice permutation condition on i. It follows 
that the map just described is a bijection between semi-standard Young tableaux 




where 




for i — 1, 2, . . . , r. 
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Figure 1: Corresponding semi-standard Young and Littlewood-Richardson 
tableaux 

with shape A and content /i and Littlewood-Richardson tableaux with shape 
r/cr and content A. Thus, computing Kostka numbers reduces to computing 
Littlewood-Richardson coefficients. 

As a corollary, we get a similar result for the computation of Kostka numbers 
based on polyhedral methods. 

Corollary 2.7. For fixed rank r, the Kostka number Kxfj, can be computed in 
polynomial time in the size of the input weights A and fi. For arbitrary rank 
one can decide in polynomial time in the size of the weights whether Kxfi = 
or not. 

3 Using the Algorithm in Practice 

Using the explicit definitions for the hive and BZ-polytopes as the sets of solu- 
tions to systems of linear inequalities and equalities, we wrote a Maple notebook 
which, when given a triple of highest weights, produces the corresponding hive 
or BZ-polytope in LattE readable input format. The notebook is available from 
|http : //m ath . ucdavis . edu/~tmcal All computations were done on a Linux 
PC with a 2 GHZ CPU and 4 Gigabytes of memory. 

From our experiments, we conclude that (1) The polyhedral method of com- 
puting tensor product multiplicities complements the method employed in LiE. 
LiE is effective for larger ranks (up to r = 10, say), but the sizes of the weights 
must be kept small. This is because LiE uses the Klimyk formula to generate the 
entire direct sum decomposition of the tensor product, after which it dispenses 
with all but the single desired term. However, computing all of the terms in 
the direct sum decomposition is not feasible when the sizes of the entries in the 
weights grow into the 100s. On the other hand, (2) Lattice point enumeration is 
often effective for very large weights (in particular, the algorithm is suitable for 
investigating the stretching properties of Section 0)). However, the rank must 
be relatively low (roughly r < 6) because lattice point enumeration complexity 
grows exponentially in the dimension of the polytope, and the dimensions of 
these polytopes grow quadratically with the rank of the Lie algebra. Together, 
the two algorithms cover a larger range of problems. 

We would also like to mention that other authors are also using lattice points 
in polytopes to compute Clebsch-Gordan coefficients, this time via the calcu- 
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lation of vector partition functions for classical root systems These 
authors report that, like us, they can compute with large sizes of weight entries. 

3.1 Experiments for type Ar 

In the tables below, we express highest weights in terms of the canonical basis 
ei, . . . , e^, so that the highest weights are partitions with r parts. 

Experiments indicate that lattice point enumeration is very efficient for com- 
puting Littlcwood-Richardson numbers when r < 5. First, wc computed over 
30 instances with randomly generated weights with leading entries larger than 
40 with our approach and with LiE. In all cases our algorithm was faster. After 
that, we did a "worst case" sampling for the table in Figure [5] comparing the 
computation times of LattE and LiE. To produce the ith row of that table, 
we selected uniformly at random 1000 triples of weights (A, /i, i^) in which the 
largest parts of A and /i were bounded above by lOi and \v\ = |A| + (this 
is a necessary condition for ^ ^ 0). Then we evaluated the corresponding 
hive polytopes with LattE. The LattE input files are created with our Maple 
program. The weight triple in the ith row is the one that LattE took the longest 
time to compute. We then computed the same tensor product multiplicity with 
LiE. The table in Figure |31 shows the running time needed when using LattE 
to compute weight triples with entries in the thousands or millions. 

When r > 6, the running time under LattE begins to blow up. Still, for 
r = 6, all examples we attempted could be computed in under 30 minutes 
using LattE, and most could be computed in under 5 minutes. For example, 
among 54 nonempty hive polytopes chosen uniformly at random among those 
in which the weights had entries less than 100, all but seven could be computed 
in under 5 minutes with LattE, and the remaining seven could all be computed 
in under 30 minutes. None of these Littlewood-Richardson coefficients could be 
computed with LiE. At r = 7, lattice point enumeration becomes less effective, 
with examples typically taking several hours or more to evaluate. 

3.2 Experiments for types B^, Cr, and Dr 

To compute Clebsch-Gordan coefficients in types Br, Cr, and Dr, we used the 
BZ-polytopes. In the tables that follow, all weights are given in the basis of 
fundamental weights for the corresponding Lie algebra. 

Our experiments followed the same process we used for Ar'. First, for each 
root system, we computed over 30 instances with randomly generated weights 
with entries larger than 40 with our approach and with LiE. In all cases our 
algorithm was faster. After that, we did a "worst case" sampling to produce the 
tables in Figure ^ comparing the computation times of LattE and LiE. As in 
Section IXn these weight triples were the ones which LattE took the longest to 
evaluate among thousands of instances generated with the following procedure: 
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A, fi, v 




LattE runtime 


LiE runtime 


(9,7,3,0,0), (9,9,3,2,0), (10,9,9,8,6) 


2 


0m00.74s 


OmOO.Ols 


(18,11,9,4,2), (20,17,9,4,0), (26,25,19,16,8) 


453 


0m03.86s 


0m00.12s 


(30,24,17,10,2), (27,23,13,8,2), (47,36,33,29,11) 


5231 


0m05.21s 


0m02.71s 


(38,27,14,4,2), (35,26,16,11,2), (58,49,29,26,13) 


16784 


0m06.33s 


0m25.31s 


(47,44,25, 12, 10), (40,34,25, 15,8), (77,68, 55,31,29) 


5449 


0m04.35s 


lm55.83s 


(60,35, 19, 12, 10), (60,54,27,25,3), (96,83,61,42,23) 


13637 


0m04.32s 


23m32.10s 


(64,30,27,17,9), (55,48,32,12,4), (84,75,66,49,24) 


49307 


0m04.63s 


45m52.61s 


(73,58,41,21,4), (77,61,46,27,1), (124,117,71,52,45) 


557744 


0m07.02s 


> 24 hours 



Figure 2: A sample eomparison of running times between LattE and LiE case of 







LattE runtime 


(935,639,283,75,48) 

(921,683,386,136,21) 

(1529,1142,743,488,225) 


1303088213330 


0m07.84s 


(6797,5843,4136,2770,707) 
(6071,5175,4035,1169,135) 
(10527,9398,8040,5803,3070) 


459072901240524338 


0m09.63s 


(859647,444276,283294,33686,24714) 
(482907,437967,280801,79229,26997) 
(1120207,699019,624861,351784,157647) 


11711220003870071391294871475 


0m08.15s 



Figure 3: Computing large weights with LattE. Case of Ar 



First, to produce line i of a table, we set an upper bound Ui for the entries of 
each weight. Then, we generated 1000 random weight triples with entry sizes no 
larger than Ui. Here are the specific values of Ui used in each of the three tables 
in Figure^ For type Br, the bounds Ui were 50, 60, 70, and 10,000, respectively. 
For type Cr, the bounds Ui were 50, 60, 80, and 10,000, respectively. Finally, 
for type Dr, the bounds Ui were 20, 30, 40, and 10,000, respectively. For each 
generated triple of weights, we produced the associated BZ-polytopes (using our 
Maple notebook) and counted their lattice points with LattE. In the table are 
those instances that were slowest in LattE. We also recorded in the table the 
time taken by LiE for the same instances. One can see the running time needed 
by LattE is hardly affected by growth in the size of the input weights, while the 
time needed by LiE grows rapidly. 

Wc found that for types Bj. and Cr, lattice point enumeration with the BZ- 
polytopes is very effective when r < 3. Each of the many thousands of examples 
we generated could be evaluated by LattE in under 10 seconds (the examples 
in Figure0]were the worst cases). When r ~ 4, the computation time begins to 
blow up, with examples typically taking half an hour or more to compute. The 
polyhedral method is also reasonably efficient for type D Lie algebras with rank 
4, the lowest rank in which they are defined. All of the examples we generated 
could be evaluated by LattE in under 5 minutes. 

4 Two Conjectures that Could Generalize the 
Saturation Theorem 

In 1998, Knutson and Tao used the hive polytopes to prove the Saturation 
Theorem. Buch has written a very clear exposition of this proof in . 

Theorem 4.1 f |19p. (Saturation) Given highest weights A, /i, and v for a 

Lie algebra of type Ar, and given an integer n > Q, the Littlewood-Richardson 
coefficient c'^ satisfies 



In hive polytopc language, the Saturation Theorem can be restated as 






where d = . The definition of hive polytopes (see Definition 12.31 above) 
impHes that HnXnfj, — '^-^A/^' Saturation Theorem is equivalent to the 

existence of a lattice point inside each hive polytope H^^. Now we would like 
to state two conjectures that generalize this theorem. 
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LattE runtime 


LiE runtime 




(46,42,38), (38,36,42), (41,36,44) 
(46,42,41), (14,58,17), (50,54,38) 
(15,60,67), (58,70,52), (57,38,63) 
(5567,2146,6241), (6932,1819,8227), (3538,4733,3648) 


354440672 
88429965 
626863031 
215676881876569849679 


0m09.58s 
0m06.38s 
0m07.14s 
0m7.07s 


lm45.27s 
3ml6.01s 
6m01.43s 

n/a 


Cz 


(25,42,22), (36,38,50), (31,33,48) 
(34,56,36), (44,51,49), (37,51,54) 
(39,64,58), (65,15,72), (70,41,44) 
(5046,5267,7266), (7091,3228,9528), (9655,7698,2728) 


87348857 
606746767 
519379044 
1578943284716032240384 


0m07.48s 
0ni08.42s 
0m07.63s 
0ni07.66s 


0ml7.21s 
2m57.27s 
8m00.35s 
n/a 




(13,20,10,14), (10,20,13,20), (5,11,15,18) 
(12,22,9,30), (28, 14,15,26), (10,24,10,26) 
(37,16,31,29), (40,18,35,41), (36,27,19,37) 

(2883,8198,3874,5423), (1901,9609,889,4288), (5284,9031,2959, 5527) 


41336415 
322610723 
18538329184 
1891293256704574356565149344 


2m46.88s 
3m04.31s 
4m29.63s 
2m06.42s 


0ml2.29s 
7m03.44s 
>60m 
n/a 



Figure 4: A sample comparison of running times between LattE and LiE 



4.1 First Conjecture 



To show that every hive polytope contains an integral point, Knutson and Tao 
actually proved that every hive polytope contains an integral vertex. Our idea 
was to take a different approach to show a generalization of this last result 
using the basic theory of triangulations of semigroups. To develop this idea, 
observe that the boundary equalities and rhombus inequalities that define a 
hive polytope may be expressed as the set of solutions to a system of matrix 
equalities and inequalities: 



h G 



»(r+l)(r+2)/2 



Bh = b{X, /i, I/), 
Rh < 



(2) 



where B and R arc integral matrices (depending on ?■), and b{X, fi, v) is a integral 
vector depending on A, fi, and v. Here we think of a hive pattern h as a column 
vector of dimension (r + l)(r + 2)/2. Note that there is some degree of choice in 
how the boundary equalities and rhombus inequalities are encoded as matrices B 
and R, respectively. However, all such encodings arc equivalent for our purposes. 

A polytope defined by such a system of equalities and inequalities may be 
homogenized by adding "slack variables" . This produces an equivalent poly- 
tope defined as the set of nonnegative solutions to a system linear equalities. 
Following this procedure, we define the homogenized hive polytope H^^ by 



HI, = i h 



B 
R I 



h = 



6(A,^, 1^) 




h>0 



(where / is the identity matrix). The equivalence between H^^ and H^^ is given 
by the linear map 

h 

-Rh ■ 



h 



Note that this linear map preserves vertices and integrality. Therefore, to prove 
the Saturation Theorem, it suffices to show that every homogenized hive poly- 
tope contains an integral vertex. Proceeding with this idea, we make the fol- 
lowing definitions. 



Definition 4.2. Fix r G Z. Define the homogenized hive matrix to be 

M = 



B 
R I 



(where B and R are as in equation Given an integral vector b with dimen- 
sion equal to the number of rows in M , define the generalized hive polytope or 
g-hive polytope iJf, by 



Hb = \h: Mh ^b, h>0 



(3) 
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Note that the homogenized hive polytopes are g-hive polytopes which satisfy 
additional restrictions to the right-hand side vector (such has the final entries 
of b being 0). 

Wc now state sonic very basic facts about vertices of polyhedra expressed 
in the form {x : Ax ~ b, x > 0}. Let a finite collection of integral points 
{ai, . . . , a„} C be given, and let A be the matrix with columns oi, . . . , a„. 
Define cone^ to be the cone in R"* generated by the point-set {ai, . . . , a„}: 

cone A = {xiai + • • • + x„a„ : xi, . . . , x„ > 0}. 

Then, for each vector b £ Z™, we have a polytope 

Pb = {x : Ax = b, X > 0}, 

and P;, ^ if and only if 6 G cone A. In other words, there is a correspondence 
between nonempty polytopes Pb, b € Z™, and the elements of the semigroup 
generated by the columns of A. The crucial property for our purposes is the 
following. 

Lemma 4.3. If b E (cone^') fl Z™ for some m x m submatrix A' of A with 
det A' = ±1, then Pb has an integral vertex. 

Proof. Suppose that b E (coneA') nZ™ for some mxm submatrix A' of A with 
det A' = ±1. 

Let the columns of A' he Ui^, . . . , a.i^ , and let J = {zi, . . . , im} be the indices 
of these columns. Then there is a vector x = (xi, . . . , a;„ )^ e K^o such that 
Ax — b and Xi — for each i ^ J. Letting x' — (x^^ , . . . , x^^) and using 
Crammer's rule to solve for x' in A'x' — b, we find that x is an integral vector. 
Thus, X is an integral lattice point in the polytope Pb. 

To see that x is in fact a vertex of Pb, Recall that the codimension (with 
respect to the ambient space) of the face containing a solution to a system of 
linear equalities and inequalities is the number of linearly independent equalities 
or inequalities satisfied with equality. Observe that a:: is a solution to the system 
of n equalities 

Ax = b 

Xi = 0, i ^ J. 

We claim that this is a linearly independent system of equalities. For suppose 
otherwise. Then the zero vector is a nontrivial linear combination of the rows 
of A and the row- vectors e^, i ^ J. But this implies that the zero vector is 
a nontrivial linear combination of the rows of A' , which is impossible because 
det A' ^ 0. 

Thus, having shown that x satisfies the n linearly independent equalities 
above, we have shown that x lies in a codimension-n face of Pb, i.e., x is a 
vertex. □ 
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We say that , . . . , o^^^ is a unimodular subset if the submatrix A' of A 
with columns a^^ , . . . , 0^,^ satisfies det A' = ±1. We say that the matrix A has a 
unimodular cover (rasp, unimodular triangulation) if the point set {ai, . . . , a„} 
has a unimodular cover (rcsp. unimodular triangulation). 

Corollary 4.4. //A /las a unimodular cover, then Pf, has an integral vertex for 
every integral b £ cone(^). 

Our conjecture is that this corollary applies to the homogenized hive matrix. 
More precisely, we conjecture the following. 

Conjecture 4.5. For each r, the homogenized hive matrix has a unimodular 
triangulation. Consequently, every g-hive polytope has an integral vertex. 

Since the hive polytopes are special cases of the g-hive polytopes, Conjecture 
14.51 generalizes the Saturation Theorem. 

Theorem 4.6. Conjecture i4-5{ l is true for r < 6. 

To compute the unimodular triangulations that provide a proof of Theo- 
rem ^21 we used the software topcom j^HI- It may be worth noting that the 
triangulations used to prove Theorem 14 . 61 were all placing triangulations. 

4.2 Second Conjecture 

For our second conjecture, we looked at general properties satisfied by the tensor 
product multiplicities for semisimple Lie Algebras of types Br, Cr, and Dr under 
the operation of stretching of multiplicities. By stretching of multiplicities, we 
refer to the function e: Z>o Z>q defined by e(n) = C^An^i- 

It follows from the definitions of the BZ-polytopes that, given any highest 
weights A, /x, ^ of a semisimple Lie algebra, C"j^ „^ = e(n) is a quasi-polynomial 
in n. Indeed, e(n) is, in polyhedral language, the Ehrhart quasi-polynomial of 
the corresponding BZ-polytope. We recall the basic theory of Ehrhart quasi- 
polynomials. Its origins can be traced to the work of Ehrhart ^Hl in the 1960's 
(see Chapter 4 of for an excellent introduction). 

Given a convex polytope P, let 



If P is a d-dimcnsional rational polytope in M'"', then the counting function 
ip{n) = ff{nP n Z*"') is a quasi-polynomial function of degree d; that is, there 
are polynomials fi(n), . . . , fN{n) of degree d s.t. 



nP^{x: {l/n)x E P} 



n = l,2,... . 
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If we put P = H^^, then the Ehrhart quasi-polynomial of P is just the 
stretched Littlewood-Richardson coefficient cj^^l.n/i- The Ehrhart quasi-polyno- 
mials of hive polytopes have been studied by several authors. Our experiments 
support the conjecture of King et al jEl that the coefficients are in fact positive. 
Since large weights can be computed with lattice point enumeration, it is possi- 
ble to produce the Ehrhart quasi-polynomials for the stretched Clebsch-Gordan 
coefficients in the other types. See Figures |SH3 for some sample examples out of 
the many hundreds generated. Our experiments motivate the following "stretch- 
ing conjecture" . 

Conjecture 4.7. (Stretching Conjecture) Given highest weights X^^^v of a Lie 
algebra of type A.^, Bj., C,., or D^., let 

{fi{n) if n = 1 modN, 
■ 
fN{n) ifn = N modN 

be the quasi-polynomial representation of the stretched Clebsch-Gordan coeffi- 
cient C^^ „^ . Then the coefficients of each polynomial fi are all nonnegative. 

The type case of this conjecture was made by King, ToUu, and Toumazet 
in ^H]. That Conjecture 14.71 implies the Saturation Theorem follows from a 
result of Derksen and Weyman showing that the Ehrhart quasi-polynomials 
of Hive polytopes are in fact just polynomials. 

We should remark that the saturation property is known not to hold in the 
root systems i?r, Cr, and Dr- A simple example in B2, due to Kapovich, Leeb, 
and Millson is given by setting A = /x = ;/=(l,0) (with respect to the 
basis of fundamental weights) . In this case we have 

^ jo if n is odd, 
II II n IS even. 

This example also demonstrates why the Stretching Conjecture is not contra- 
dicted by the failure of the saturation property in the root systems B^, Cr, or 
Dj.- Since the stretched multiplicities arc not necessarily polynomials in these 
cases, it is possible for them to evaluate to zero for some nonnegative integer 
while still having all nonnegative coefficients. 

Acknowledgements: We are thankful to Arkady Berenstein, Charles Co- 
chet, Misha Kapovich, Allen Knutson, Peter Littclmann, Ezra Miller, Hariharan 
Narayanan, Jorg Rambau, Francisco Santos, Anna Schilling, Etienne Rassart, 
Michele Vergne, and Andrei Zelevinsky for useful comments. 
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A, /i, V 




LattE runtime 


(0, 15, 5) 
(12, 15, 3) 
(6, 15, 6) 


\ 


\ 68339 ^5 _^ 407513 „4 + 13405 ^3 + 9499 „2 + 107 „ _^ ^ ^^^^ 

68339 „5 1 407513 „4 , 13405 „3 , 16355 „2 , 659 „ , 75 „ „J j 
L 64 ' 384 " ' 32 " ' 192 64 128' " 


7m08.57s 


(4, 8, 11) 
(3, 15, 10) 
(10, 1, 3) 


\ 


^ + 3 ?i + 1 , n even 
+ 3?i. + 3/4, n odd 


0m01.64s 


(8, 1, 3) 
(11, 13, 3) 
(8, 6, 14) 


\ 


' 121 6 , 1129 5 , 6809 „4 , 163 ^3 , 2771 2 , 191 , 
576 ^ 640 ^ 1152" ^ 16 ^ 288 ^ 40 " ^ -^' ^-vcii 
121 6 I 1129 5 , 6809 4 , 1933 3 , 659 2 , 8003 „ _i 93 j j 

.576 " ^ 640 ^ 1152 " ^ 192 " ^ 72 " ^ 1920 ^ 128' " 


0ml0.26s 


(8, 9, 14) 
(8, 4, 5) 
(1, 5, 15) 


\ 


f '92' + '64o' + "'f + + If + "'o' " + 1' « even 

4117 6 1 50369 „5 , 14829 „4 , 5599 „3 , 3451 „2 , 5001 „ 1 97 j j 
( 192 " 640 " ' 128 64 " 96 " 640 " ^ 128 ' " "^^^ 


0ml3.29s 


(10, 5, 6) 
(5, 4, 10) 
(0, 7, 12) 


\ 


f 669989 „5 : 286355 „4 , 10803 „3 , 7993 „ 2 , 1427 „ , -1 

960 ^ 384 ' 32 " ^ 96 120 ""1"^' t^veii 
669989 5 , 286355 4 , 10803 3 , 15509 2 , 10081 „ 1 65 j j 
L 960 ' 384 " ' 32 " ' 192 " ' 960 "^128' 


2m52.39s 



Figure 5: Stretched Clebsch-Gordan coefficients for B-^. 



A, /X, V 




LattE runtime 


(1,13,6) 

(14,15,5) 

(5,11,7) 


\ 


( 5937739 „6 , 87023 „5 , 936097 4 , 27961 „3 , 85397 2 , 883 „ , i „ 

5760 ' 40 " ' 576 " ' 48 ' 720 60 " ^^"^^^ 
5937739 ^6 , 87023 ^5 , 936097 ^4 , 27961 ^3 , 657931 ^2 , 3097 ^ i / a :„ ^AA 
[ 5760 ' 40 ' 576 " ' 48 ' 5760 240 + •J/ ft uuu 


21m20.59s 


(4,15,14) 

(12,12,10) 

(4,9,8) 


\ 
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[ 2880 " ' 960 "' ' 1152 "' ' 32 " ' 1440 " ' 960 " 128' "■ 


17m05.74s 


(9,0,8) 

(8,12,9) 

(7,7,3) 


1/30 + 3/8 + yI n3 + f n2 + „ + 1 


OmOO.eis 


(10,2,7) 
(8,10,1) 
(7,5,5) 


\ 


f 596153 „6 1 53425 „5 , 502621 4 , 5577 3 , 11941 2 , 149 „ , i „ „,„^ 
1152 "' ' 48 ' 576 ' 16 ' 144 12 -^' "■ Lveii 
596153 „6 , 53425 „5 , 502621 4 , 5577 3 , 94097 „2 , 131 „ , 23 „ „ j j 
( 1152 i 48 " ' 576 ' 16 " ' 1152 '* 12 "'"'"32' " 


19m24.55s 


(10,10,15) 

(11,3,15) 

(10,7,15) 


\ 


f 6084163 ^6 + 507527 „5 _^ 1185853 „4 + 59995 „3 + 43039 ^2 357 ^ + 1^ „ even 

6084163 „6 1 507527 „5 , 1185853 „4 , 59995 „3 , 144751 „2 , 883 „ , 25 „ _ j j 
[ 320 " ' 30 ' 192 ' 48 " ' 960 " "*" 80 " """ 64 ' " 


16m05.08s 



Figure 6: Stretched Clcbsch-Gordan coefficients for C3. 
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LattE runtime 
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\ 
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(1, 12, 0, 3) 
(0, 5, 3, 4) 


5 71^ + 4 ?i + 1 
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(12, 2, 5, 13) 
(15, 6, 10, 11) 
(2, 0, 12, 13) 




f 455263 „7 , 447281 „6 , 198433 „5 ; 971011 „4 , 108787 3 , 28969 „ 2 , 12631 „ , i 
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Figure 7: Stretched Clcbsch-Gordan coefficients for D^. 
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